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Abstract. The analytic formulas for structure constants of su(2S + 1) algebra in terms of 
3jm and 6j symbols of su(2) have been derived for the decomplexification of the Liouville- 
von Neumann equation. 



1. Introduction 

Let {Cf , Cf , C%} be a base of su(2S+l) algebra, where S = 1/2, 1, 3/2, . 
is the spin quantum number, n = (2S + l) 2 — 1. We have according to [T] 



Cfcf = - d E5 i3 + zf jh Cl Tr Cf = 0, TkCfCf = a%, (1) 



hence 



z ijk ~ 9ijk ^ e ijki (2) 

- i[Cf, Cf] = 2eg- fc Cf , {Cf , Cf} = C ~Eki + 2^^, (3) 
e? jk = ±-Tr [Cf,Cf]C%, (4) 

4 fc = lTV{Cf,C/}Cf, (5) 

where d = 2S + 1, i? is the unit matrix in dimension d x d, c is a constant, 
Tr is a symbol for trace. It is easy to see that the structure constants efj k 
and gfj k are completely antisymmetric and symmetric in the displacement of 
any pair of indices. 

2. Hermitian basis 

In order to calculate the structure constants we have to choose the basis. 
The basis is based on linear combinations of irreducible tensor operators. 



Exemplary DSID style 
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The matrix representations of irreducible tensor operators T kq [2] can be 
calculated using the Wigner 3jm symbols: 

s 

T k s q = ^/(2S + l)(2k + l) (-l) S - m (fm k q i')\S,m><S,m% (6) 

m,m'=—S 

where < k < 25, and — k < q < k in steps of 1. The normalization is 
such that Tq = E. It is known that the Cartesian product operators S x ,S y , 
and S z for spin 5 = \ are Hermitian and can be calculated from irreducible 
tensor operators [3] 

i 11 i 1 /o l\ 
S^=^= 2 (T lU -T^ = - ( V (7) 



* :(rL + ri) = i(? (8) 



Allard and Hard [4j have formed linear combinations of the irreducible tensor 
operators not only for single-quantum coherences, but also for all coherences 
according to 



C, 



S k , gx = t/^^Ct-, + (-1)*Z&), q + 0, (10) 



C S km = iJ ^^(T fc % - (-iyT k s q ), q + 0, (11) 



CL = \\ S -^pt T S q = ,k>l, (12) 



C^ = j^E, (13) 

where 1 < k < 25", and 1 < q < k in Cg ?a; , C| gy and 1 < k < 25, in C| z in 
steps of 1. The matrices (I10|ll|12j) are traceless and their number is equal 
to (25 + l) 2 — 1. Using the well-known relations for the irreducible tensor 
operators 

(T k s q ) + = (-l)«T fc ,_ g! (14) 

we can see that matrices f)10|l 1 112[) are Hermitian. Using the formula from 


TrTg jq Tg l>ql = (-1)9(25 + 1)6^6^ (15) 
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it is easy to show that the basis is normalized so that S x = Cf x , S y = Cf y , 
S z = Cf z , irrespective of the spin quantum number 5, i.e. 



(C r , C s ) — Tr C r C s — 5 r 



5(5 + 1)(25 + 1) 



(16) 



The set (|10til3j) is complete. The matrices Ck, z are diagonal 

[Ct,C§J=0. (17) 

There also exist the other useful bases [5]. Prom the physical point of view, 
for some applications the basis [3] is preferred. 



3. Analytic formulas for structure constants 

There are 27 combinations in threes including the repetitions: XX' X' 

^k'.a'xi 1 — ^k",q"yi ^ ' 



XX'Y", XX' Z". . . , where X 



C^n z and so on. The use of the symmetrical properties of the Wigner 3jm 
symbols and the formula 2.4.(23) [2] 

TrTS q Tg w Tg, w , = (-l) 25 + fc + fc '+ fc "(25 + 1)§ 
[(2* + l)(2k' + l)(2k" + !)]*{* k s I"} (J $f) , 



(18) 



where { s s g } is the 6j symbol, allows us after substitution of (|10lllll2p 
in (|4|5h , to calculate all structure constants of su(2S + 1) algebra. Let us 
introduce the function 



k k' k"\ 
S S S I- 



(~l) 2S I 

F(k, k , k", S) = ^-^JS(S + 1)(25 + l)(2k + l)(2k' + l)(2k" + 1) { 
V3 

(19) 

All antisymmetric structure constants are zero for K = k + k' + k" even 
and nonvanishing antisymmetric structure constants in terms of 3jm and 6j 
symbols have the explicit form are presented by formulas (|20|21l22p for K 
odd: 



P s 

e XX'Y" 



F r 

72 



-iy (;_*; - fe ;») + (-¥ (*, 5:4") + h) 9 " ft 



(20) 



F 



(k k' k"\ 

{-q q> q" ) 



+ 



(21) 



e S _ w( iw (k k' k"\ 

e XY'Z" - ~* [q -q> ) 



(22) 
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All symmetric structure constants are zero for K = k + k! + k" odd and 
nonvanishing symmetric structure constants in terms of 3jm and 6j symbols 
have the explicit form are presented by formulas (|23|24I25|) for K even: 

(23) 



9xx'X" 



r 

777 



dXY'Y" 



F r 

7i 



A; fc' fc'" 



(24) 



«W = = ^(-l) 9 (} 4' o") , fffz'Z'' = ^(-l) 9 (8 o'o") • (25) 

We have in X, Y 1 < k, k' , k" < 25, 1 < q < k, 1 < q' < k', 1 < q" < k" and 
in Z 1 < k, k', k" < 25 in steps of 1. 

The direct calculation confirms that the structure constants efj k and gf^ k are 
completely antisymmetric and symmetric in the displacement of any pair of 

operators. In other words it is e-xx'Y" = ~ e XY"X'i 9xx'Z" = 9yy'Z" an d so 
on. 



4. Decomplexification of the Liouville-von Neumann equation 

The structure constants of su(2S + 1) algebra have wide physical appli- 
cations. Let us consider the Liouville-von Neumann equation for the density 
matrix p, describing the dynamics of a system. It has the form 

idtp = [H,p], p(t = 0) = po, p + = p, Trp = 1, (26) 

where H is the Hamiltonian of the system. 

The state of a cf-dimensional quantum system (qudit) is usually described 

by the dxd density matrix p. Let us present the solution of the equation (|26|) 

for one qudit as p = -, R a C^, R = 1, H = ^hpCg 1 . We 

(25i+l)y Sl(s 3 1+1) 

multiply (|26|) by C^ 1 and execute the operation Tr , where by definition 

TrpC* = ^/^f^R a , Ct E {C^.,Cf w; C*,c£}. Hereinafter, the 
summation is made over the repeating Greek indices on the set (jlOM13p . and 
over Latin indices on the set (|10lllll2p and then on k, and then on q. The 
Liouville-von Neumann equation takes on the real form in terms of the func- 
tions Rj as a closed system of differential equations for the set of initial 
conditions [6]: 

dtRi = ef^hiRj. (27) 
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The length of the generalized Bloch vector b Sl is conserved under unitary 
evolution 

b S > = y/^. (28) 

For two different coupled qudits we have H = ^hapC^ 1 ® C^ 2 , 

9 = Wi + ms 2 + vTWi + + tj^ 1 Cs 2 ' (29) 

i^oo = 1) where (g) is the symbol of direct product. 

The dynamic equation (|26p for two different qudits takes on the real form 
in terms of the functions R m Q, Ro m , R mn as a closed system of differential 
equations (|3QI31l32p for the set of initial conditions: 



d t R 



mO 



'5 2 (5 2 + l) 



(30) 



d t R 



0/71 



/ 5i(5i + l) „ Sa 

"pirn 



<"'. -—(hopRoi + hipRn) 



(31) 



t^mn 



e Sl 

pirn 
"pin 



S 2 {S 2 +l) 
3 



(h pn Rio + h p oRi n ) + g^f n h pr R. 



+ 



SiQSi+i) 

3 



(hmpRoi + hopRrai) + gfi m h rp Ri 



where by definition 



(32) 



1 



TrpCf 1 ® Cf 2 = ^5i(5i + l)S 2 (S 2 + l)i? a/3 , 



(33) 



and Cf 2 G {^.C^C^Cg} . The functions R m o,Rom describe the 
individual qudits and the functions R mn define their correlations. 

The length of the generalized Bloch vector b Sl 2 is conserved under uni- 
tary evolution 

foSlS2 _ 



7?2 i r>2 i p2 
n m0 "i" n 0m ^mn- 



(34) 



The set of equations for 3 qubits has been obtained in [7J. The dynamic 
equations for 3 different qudits will be presented elsewhere. 



5. Conclusion 

It is not necessary for the basis to be Hermitian since the results of calcu- 
lations are independent of the choice of base, but there is the main advantage 
with the Hermitian basis. It is that the Liouville-von Neuman equation not 
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involve any complex numbers and can be solved using real algebra. This is 
not true for non-Hermitian bases. Real algebra makes numerical calculations 
faster and simplifies the interpretation of the equation system. 
This basis forms a natural basis for calculations on coupled spin systems 
(multipartite systems) [4j because all the single-spin operators are part of 
the complete basis when the unit operator is part of the single-spin basis. 
The convolutions of structure constants [8] can give rise to the additional 
relations between 3jm and 6j symbols. 

The relationship between the basis [4] and the Gell-Mann basis for spin 1 has 
been presented in the appendix. 



Appendix 



The matrix representation of the complete set of the Hermitian operators 
for spin 1 (for qutrit) has the form 




r 1 



C 1 



C 1 



These matrices are traceless TrC a = and orthogonal Tr C a Cb = 25 a b, 
1 < a, b < 8. The relationship between the basis C a [3] and the Gell-Mann 
basis A a is the following: 

d = l/V2(Ai + A 6 ),C 2 = 1/V2(A 2 + X 7 ),C 3 = 1/2A 3 + s/E/2X 8 ,C 4 = 
A 5 , C 5 = 1/V2(A 2 - A 7 ), C 6 = V3/2A 3 - 1/2A 8 , C 7 = l/ v / 2(Ai - A 6 ), C 8 = A 4 . 
The antisymmetric (symmetric) structure constants e a b c (g a bc) are corre- 
spondingly equal to: ei 2 3 = ei 58 = e 254 = e 278 = e 375 = e 47 i = 1/2, ei 56 = 

e672 = V / 3/2, e 348 = —1; g 33 6 = <? 44 6 = —5666 = 5688 = l/\/3, 5556 = 5116 = 
5226 = 5677 = — 1/2a/3, 5235 = 5118 = 5558 = 5124 = 5137 = ~ 5228 = ~ 5778 = 

-5475 = 1/2. Hence, we have C a C b = 2/3<5 ab + g a b c C c + ie a b c C c - 
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